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We provide exact expressions for the electrostatic energy of uniformly-charged prolate and oblate spheroidal 
shells. We find that uniformly-charged prolate spheroids of eccentricity greater than 0.9 have lower Coulomb 
energy than a sphere of the same area. For the volume-constrained case, we find that a sphere has the high¬ 
est Coulomb energy among all spheroidal shells. Further, we derive the change in the Coulomb energy of a 
uniformly-charged shell due to small, area-conserving perturbations on the spherical shape. Our perturbation 
calculations show that buckling-type deformations on a sphere can lower the Coulomb energy. Finally, we con¬ 
sider the possibility of counterion condensation on the spheroidal shell surface. We employ a Manning-Oosawa 
two-state model approximation to evaluate the renormalized charge and analyze the behavior of the equilibrium 
free energy as a function of the shell’s aspect ratio for both area-constrained and volume-constrained cases. 
Counterion condensation is seen to favor the formation of spheroidal structures over a sphere of equal area for 
high values of shell volume fractions. 


I. INTRODUCTION 

Shapes of physical systems as diverse as blood cell mem¬ 
branes, colloidal particles, nanowires, and galaxies are often 
considered as spheroidal with varying degrees of eccentric¬ 
ity. Many charged structures such as colloids or emulsions are 
usually modeled as spheroidal shells with a uniform surface 
charge density. To our best knowledge, the expression for the 
Coulomb energy of a uniformly-charged spheroidal shell is 
not available in the literature. In this article, we provide the 
needed result. Further, we derive a general expression for the 
change in the Coulomb energy of a uniformly-charged shell 
due to small, area-conserving perturbations on the spherical 
shape. Using the result, we explore the existence of defor¬ 
mations that can lower the electrostatic energy relative to the 
unperturbed charged sphere. 

We note that a closed-form expression for the electrostatic 
potential energy of a solid homogeneously-charged spheroid 
can been obtained [1]. Further, the electrostatic energy of a 
conducting spheroidal shell as a function of the aspect ratio is 
available elsewhere [2]. Calculations of the electrostatic po¬ 
tential for a system of point charges inside dielectric spheroids 
[3, 4] have been performed as well. On the other hand, it 
is useful to note the work in the context of solving Poisson- 
Boltzmann equation in spheroidal geometry [5, 6]. In this pa¬ 
per, we provide a comprehensive study of the homogeneously- 
charged spheroidal shell system which has been missing in the 
literature. 

Primary motivations behind our calculations stem from the 
study carried out in Ref. 7, where equilibrium shapes of 
charged, soft shells constrained to maintain a fixed volume 
were analyzed using molecular dynamics simulations. Some 
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of the results derived in the present article were employed 
to verify the oblate-shaped shell structures found in Ref. 7 
and calculate the effects of ion condensation on the equilib¬ 
rium shape of these structures. The supporting information 
associated with Ref. 7 also contained a brief derivation of the 
electrostatic energy of a uniformly-charged oblate spheroidal 
shell. In this paper, we derive the general expression for the 
Coulomb energy of prolate spheroidal shells. For the sake of 
completeness we also include the derivation of the Coulomb 
energy of oblate shells showing details that were omitted in 
Ref. 7. We analyze the variation of the Coulomb energy 
of spheroidal shells, subject to the constraints of fixed area 
or volume, as the aspect ratio of the shell is changed. We 
also examine the effects of ion condensation, computed via a 
Manning-Oosawa two-state model analysis [8, 9], on the vari¬ 
ation of the equilibrium free energy of the shell-counterion 
system. Finally, we note that it is straightforward to augment 
the energy expressions obtained here to reveal the gravita¬ 
tional potential energy of a uniformly-dense spheroidal sur¬ 
face which is often used as a model to study galaxies [10]. 

The key findings of this paper are: i) A homogeneous pro¬ 
late (cigar-shaped) spheroidal shell with eccentricity greater 
than ^ 0.9 has a lower electrostatic energy than a spherical 
shell of the same area. The lowest-energy shape of the shell, 
constrained to maintain its area, is a very long and thin pro¬ 
late spheroid whose energy approaches zero as its major-axis 
length is stretched to infinity, ii) For shells that are constrained 
to maintain their volume, the spherical shape has the maxi¬ 
mum Coulomb energy. An infinitely long and thin wire-like 
shape and a thin, flat disc of infinite area, are the degenerate 
lowest-energy shapes with vanishing energy, iii) Perturbation 
calculations show that, for the case of fixed area constraint, the 
Coulomb energy of a uniformly-charged sphere can be low¬ 
ered by a buckling-type deformation, iv) Results from the 
two-state model approximation of the shell-counterion sys¬ 
tem show that counterion condensation favors the formation 
of spheroidal structures over a sphere of equal area for high 
values of shell volume fractions. 
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The paper is organized as follows. In Sec. II, we provide 
the expression for the Coulomb energy of a uniformly-charged 
spheroidal shell, discuss the important limiting cases, and spe¬ 
cialize the expression for the case of constant area and con¬ 
stant volume constraints. Sec. Ill shows the comparison be¬ 
tween the energy of a sphere and nearly-spherical structures 
formed by a small, generic perturbation around the spherical 
shape. In Sec. IV, we discuss the effects of charge renormal¬ 
ization on the energies obtained for the spheroidal shell sys¬ 
tem and Sec. V is the conclusion. Appendices A and B present 
the derivation of the Coulomb energy of uniformly-charged 
prolate and oblate spheroidal shells respectively, and in Ap¬ 
pendix C we derive the electrostatic energy of a uniformly- 
charged circular disc. 
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FIG. 1. Cross-section of a spheroid normal to the j/-axis with c and 
a as the semi-axis lengths and A = c/a being the spheroid’s aspect 
ratio, (a) Prolate spheroid with A > 1. (b) Oblate spheroid with 
A < 1. 


II. ELECTROSTATIC ENERGY OE 
UNIEORMLY-CHARGED SPHEROIDAL SHELLS 

Consider a spheroidal shell with charge Q distributed uni¬ 
formly over its surface such that the charge density is given by 
a = Q/A, where A is the area of the spheroid. A spheroid is 
an ellipsoid two of whose semi-principal axes are equal. As¬ 
suming that the equal lengths correspond to the dimensions 
along the x and y axes, such that the cross section normal to 
the z-axis is a circle, we can describe the spheroidal shell via 
the equation; 



where a and c are the semi-principal axes, and r = y/+ y'^ 
is the distance between the point on the surface of the spheroid 
and the z-axis. A prolate spheroid is a spheroid where c > a, 
whereas an oblate spheroid corresponds to the c < a condition 
(see Fig. 1). Clearly when a = c, the spheroid reduces to a 
sphere. 

It is convenient to characterize the spheroid by defining the 
aspect ratio A defined as 

A=-. (2) 

a 

Values of A < 1 correspond to oblate spheroid whereas a pro¬ 
late spheroid is associated with A > 1. A —?► 0 corresponds to 
a circular disc, A = 1 is a sphere, and A —oo limit produces 
an infinitely long and thin rod-like spheroid. It is also useful 
to introduce the eccentricity Cp of a prolate spheroid defined 
as: 


leads to a circular disc and ^ 1 corresponds to a very long 
and thin rod like shape. We will invoke these limits at several 
places in what follows. 

In Appendix A, we derive the electrostatic energy of a 
uniformly-charged prolate spheroidal shell. Our calculations, 
which employ the standard method of separation of variables 
[11, 12], lead to the following result: 

t/(ep, c, tr) = dTT^cr^c^-x 

ep 

^ ?!^P„(l/ep)g„(l/ep)(i^„(ep))^ 

nGeven 

( 5 ) 

where n is an even integer, and Qn are Legendre functions 
of the first and second kind respectively, and Hn{e^) is the 
integral 


Hn{ep) = j — e'^cos'^v Pn{cosv)smv dv. (6) 

It is useful to express the result in terms of the total charge Q 
rather than tr. Using the fact that the area of a prolate spheroid 
is 

.!4(ep,c) = 27rcyi-e2T(ep), (7) 

where 




Similarly, we have the eccentricity Cq for an oblate spheroid: 



Note that either eccentricities lie between 0 and 1. When 
Co, Cp —> 0, the spheroid reduces to a sphere. The limit Co ^ 1 


we obtain the electrostatic energy of a homogeneously- 
charged prolate spheroidal shell to be: 


u,= 


2c epr(ep)2 


^ (2n+l)P„ 

nGeven 



Qn 





( 9 ) 


where C/p = f/p(ep, c, Q). Eq. (9) provides the first key re¬ 
sult of this paper. We point out that the above expression 
for the Coulomb energy is obtained assuming that the pro¬ 
late spheroidal shell is in vacuum. If the medium surrounding 
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the shell is polarizable and uniform, the above result for the 
energy must be scaled down by the dielectric constant of the 
medium. 

Appendix B provides the derivation for the Coulomb en¬ 
ergy of the oblate spheroidal shell. We note that a brief ac¬ 
count of this derivation appears in the supplementary infor¬ 
mation of Ref. 7. We find the electrostatic energy to be: 


U{eo,a,a) = -x 


Go 


E 




Co 


Qn I ^ 




Go 


/n(Go)^ 


( 10 ) 


where n is an even integer and /^(go) is the integral 

In{eo) = / \Jl — e^sin^v Pnicosv)smv dv. (11) 

Jo ^ 


Once again, we express below the result in terms of the total 
charge Q. Using the fact that the area of an oblate spheroid is 

£/o{eo,a) = 27ra^5'(eo), (12) 

where 

5'(eo) = 1-b ( — - Go ) tanh'^Co, (13) 


we obtain the electrostatic energy of the uniformly-charged 
oblate spheroidal shell to be: 


C/o = 


Q'^i 


2aeoS{eoY 




Go 



(14) 


where Uq = {7o(Go, a, Q). Eq. (14) provides the second key 
result of this paper. In several physical situations, geomet¬ 
ric constraints such as the constraint of fixed area or vol¬ 
ume are naturally present and it is of interest to find the shell 
shape that minimizes the Coulomb energy when only area¬ 
preserving or volume-conserving deformations are allowed. 
Using the above expressions for the electrostatic energy, we 
analyze the Coulomb energy of a shell that is subjected to 
these constraints and present the results in sections IIB and 
IIC. Before that we take a quick look at the limiting cases of 
the energy expressions found in Eqs. (9) and (14). 


A. Limiting cases 

We recall that both the prolate and oblate eccentricities lie 
between 0 and 1. Eirst we let CpjCo approach zero, which 
corresponds to a spherical shell, and find 

t/p(Gp ^ 0,a, Q) = C/o(eo ^ 0, a, Q) = ^. (15) 


We recover the well-known result for the energy of a 
uniformly-charged spherical shell. 

Taking the limit Co —1 of the oblate energy expression in 
Eq. (14) gives 

Uo(eo ^ 1,0, Q;n < 6) = 0.84872^, (16) 

where, because of the rapid convergence of the sum in 
Eq. (14), we have included terms up to n = 6 in obtaining the 
above result. The limit Co 1 corresponds to the shape of 
a circular disc. Unlike the spherical case, an exact expression 
for the energy of a uniformly-charged disc to our best knowl¬ 
edge has not been reported in the literature. It is possible to 
arrive at this energy starting from the electrostatic potential 
on the surface of the disc derived in Ref. [13]. We show the 
derivation of the energy in Appendix C. The result is 

t/disc = —, (17) 

dTT a 

where a is the radius of the disc and Q = naJa is the to¬ 
tal charge. We can now compare the disc energy obtained in 
Eq. (16) with the above exact result and find the deviation to 
be ~ 0.01%. Clearly, the two energies are in very good agree¬ 
ment. 

The prolate energy expression, Eq. (9), in the limit of pro¬ 
late eccentricity approaching unity gives (7p(ep —>■ 1, a, Q) = 
oo, that is, the energy diverges. In this limit, the prolate 
spheroid is transformed into a thin, long rod-like shape where 
the width of the rod is shrunk at the same time as the length 
of the rod is stretched out. The divergence of the energy 
arises because, as the width narrows, the distance between the 
charges on the surface shrinks faster in comparison with the 
charges growing apart due to the extension in the length. As 
we will find out in the next section, when we impose con¬ 
straints of constant shell area or volume, the prolate energy no 
longer diverges when the aforesaid limit is taken. 


B. Charged spheroidal shells of equal area 


We begin with the results for the application of area con¬ 
straint. We will assume the reference shell shape in the analy¬ 
sis to be the sphere. The area of a prolate spheroid is given by 
Eq. (7) which we rewrite below: 

i/p(Gp,c) = Sttc^y^I - el , (ig) 

This equation suggests that if is fixed, c and Cp are coupled. 
Assuming that the area is constrained to that of a sphere of 
radius R, c and Cp become related via the equation: 


c(ep, R) = R 


\ 


1 - 
^ Gp 



(19) 


Eliminating c from Eq. (9) using the above equation, we arrive 
at the expression for the electrostatic energy of a uniformly- 
charged prolate spheroidal shell constrained to a fixed area of 
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AttB? as a function of the eccentricity Cpi 




Q- 


1 - - Isin ^Cp) 


2R^epT[epY 


X (20) 


(2n + 1) P„ (l/cp) Qn (1/ep) i?„(ep)^. 

n^even 


We proceed similarly with the case of an oblate spheroidal 
shell. The area of an oblate spheroid is given by Eq. (12) 
which we rewrite below: 


j 2 ^(eo, a) = 27ra^ 




tanh ^ 



( 21 ) 


where tanh“^ denotes the inverse hyperbolic tangent function. 
We note that in the limit Cq —S' 1, the oblate shell reduces to a 
structure resembling a circular disc having two faces with total 
area 27ra^. If s^o is fixed, Co and a are coupled and recalling 
that the area is constrained to the value AtpE?, we obtain the 
relation: 


i(eo, R) = R 


\ 1 + (^ - eo) tanh- 


( 22 ) 


Eliminating a from Eq. (14) using the above equation, we 
arrive at the expression for the electrostatic energy of a 
uniformly-charged oblate spheroidal shell subject to the con¬ 
straint of fixed area: 


Uo,a{^Qj R-/ Q) - 


Q'^i^l+ - Co) tanh 


2i?V2eoS'(eo)2 


(23) 


y] (2n + 1) 

n^even 



Qn 



It is useful to express the energies in Eqs. (20) and (23) as 
a function of the aspect ratio A defined in Eq. (2). Noting that 
Bp = v/A^ — 1/A and Bo = v/l — A^, we arrive at the result: 


C/a (A) 


U,^a{VT^,R,Q) 0<A<1 

Up,a{V>^^/X:R,Q) A>1, 


where we have suppressed the dependence of Ua on other 
variables for brevity. The values of A > 1 correspond to pro¬ 
late spheroids and 0 < A < 1 region corresponds to oblate 
spheroids. Equation (24) provides the Coulomb energy of 
uniformly-charged spheroidal shells, all having the same area, 
for values of the aspect ratio A ranging from 0 to oo. 

We now analyze the variation of Ua{X) as A is changed. 
We set Q = i? = 1 for simplicity. We first confirm that 
C/a (A —1) = 0.5 as should be the case for a uniformly- 
charged spherical shell (A = 1 ^ a = c) for the aforesaid pa¬ 
rameters. Eurther, it is easy to check that C/a (A —>■ 0) ~ 0.6 
which is equivalent to the Coulomb energy of an infinitely 
thin, uniformly-charged circular disc of unit radius. Taking 
the opposite limit, we find C/a (A — oo) = 0, which suggests 
that the Coulomb energy of a very thin and long wire-like 



FIG. 2. Coulomb energy of charged spheroidal shells of equal area as 
a function of their aspect ratio A. The region 0 < A < 1 corresponds 
to oblate shells and values of A > 1 represent prolate shells. Solid 
red line corresponds to the homogeneously-charged spheroidal shell 
and dashed green line is the result for the conducting spheroidal shell. 
The dotted blue line is energy of a spherical shell (reference) which is 
0.5 for either cases (homogeneous or conducting). The above results 
are for Q = 1, R = 1. See text for the meaning of symbols. 


shape (of finite area) vanishes. We graph the function C/a (A) 
as a function of A in Eig. 2. We choose the number of terms 
appearing in the series expansion in the energy expressions 
derived in Eqs. (20) and (23) to be n = 6 as the series con¬ 
verges rapidly. As is evident from Eig. 2, we observe that the 
sphere shape is a local minimum. However, as the aspect ratio 
is increased to values beyond ~ 2.286, which corresponds to a 
prolate eccentricity of ^ 0.9, the Coulomb energy is lowered 
below that of the sphere. The energy continues to decrease as 
A increases further and we find that the shape that corresponds 
to the lowest Coulomb energy is the very thin and long pro¬ 
late spheroidal shell (of area 47r), the minimum energy being 
0. Eurther, the energy of the oblate spheroidal shell increases 
upon increasing its eccentricity (or lowering the aspect ratio 
A), with the thin circular disc corresponding to the shape of 
maximum energy. 

It is instructive to compare the results for the case of 
uniformly-charged spheroidal shells with that of conducting 
spheroidal shells. Exact expressions for the latter have been 
obtained elsewhere [2]. Using these exact results, the particu¬ 
lar expression for the case where the area is held fixed can be 
easily derived and we summarize the final results below: 






(25) 

In the above equation, the functions a and c are given by 
Eqs. (22) and (19) respectively. Similar to the above analysis 
for homogeneously charged shells, we can evaluate the varia¬ 
tion of ^A as a function of A and we obtain the green dashed 
line in Eig. 2. We find that this line is always below the red 
solid line. This implies that allowing the surface charges to 
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move freely, as is the case with the conducting shell, lowers 
the Coulomb energy. In addition, it is important to note that 
all the prolate conducting shells have a lower energy than a 
spherical conducting shell. In sharp contrast, for uniformly- 
charged shells, the spherical shape is a clear local minimum 
as evidenced by the red solid curve in the inset of Fig. 2. 


C. Charged spheroidal shells of equal volume 

We now analyze the Coulomb energy of spheroidal shells 
that are subjected to the volume constraint, that is all the shells 
have the same volume. This analysis is very similar to the 
one presented in the last subsection and so we will keep the 
following discussion brief. The volume of a prolate spheroid 
is 


f^p(ep,c) = -TTC^ (1 - el) . 


(26) 


This equation suggests that if flp is fixed to (4/3)7ri?^, then c 
and Cp are related via the equation: 


c(ep,i?) = 


R 


(1 - 


1/3- 


(27) 


Eliminating c from Eq. (9) using the above equation, we arrive 
at the expression for the electrostatic energy of a uniformly- 
charged prolate spheroidal shell constrained to a fixed volume 
of (4/3)7ri?^ as a function of the eccentricity Cp: 


6^p,y(t'P; 7?, Q) 


2i?epr(ep)2 


(28) 


^ (2n + l)P„(l/ep)g„(l/ep)i7„(ep)2. 

nGeven 


We proceed similarly with the case of an oblate spheroidal 
shell. The volume of an oblate spheroid is given by 


f2o(eo,a) = - e^. (29) 

If rio is fixed, Co and a are coupled and recalling that the vol¬ 
ume is constrained to the value (4/3)7ri?^, we obtain the rela¬ 
tion: 


a{eo.R)= ^1/6 - (30) 

(1 - e?) 

Eliminating a from Eq. (14) using the above equation, we 
arrive at the expression for the electrostatic energy of a 
uniformly-charged oblate spheroidal shell subject to the con¬ 
straint of fixed volume: 


77o,v^(co; 77, g) 


2ReoS{eoY 


^ (2n + 1) Pn j Q 

nGeven \ ° / 



In{eof. 

( 31 ) 



FIG. 3. Coulomb energy of charged spheroidal shells of equal vol¬ 
ume as a function of their aspect ratio A. The region 0 < A < 1 
corresponds to oblate shells and values of A > 1 represent prolate 
shells. Solid red line corresponds to the homogeneously-charged 
spheroidal shell and dashed green line is the result for the conducting 
spheroidal shell. The dotted blue line is energy of a spherical shell 
(reference) which is 0.5 for either cases (homogeneous or conduct¬ 
ing). The above result is for Q = 1, i? = 1. See text for the meaning 
of symbols. 


As in the case of the area constraint, it is useful to express 
the energies in Eqs. (28) and (31) as a function of the aspect 
ratio A. Carrying out the transformation from eccentricities 
Cp, Co to A, we arrive at the result: 


7/y(A) 


C7o,y(\7l — R,Q) 0 < A < 1 
C7p.y('\/A^ — 1/A, i?, g) A > 1, 


where we have suppressed the dependence of Uy on other 
variables for brevity. Equation (32) provides the Coulomb 
energy of uniformly-charged spheroidal shells, all having the 
same volume, for values of A ranging from 0 to oo. 

In Eig. 3, we plot the change in Uy as A is varied for the pa¬ 
rameters g = 1 and R = 1. We find that Uy{X —>■ 1) = 0.5 as 
expected for the homogeneously-charged sphere. We notice 
that the sphere represents the shape of maximum Coulomb en¬ 
ergy which is in sharp contrast with the above analyzed case of 
the area-constrained charged shells. We observe in Eig. 3 that 
every other energy value is degenerate, in the sense that there 
are two distinct shapes that correspond to the same energy, one 
of the shapes being an oblate and the other a prolate. We find 
that the minimum energy for the volume-constrained charged 
spheroidal shell is 0. The two shapes that correspond to this 
value are a thin charged disc (of infinite area), and a long and 
thin prolate spheroidal shape (of infinite area). Again, this 
is in sharp contrast with the area-constrained charged shell 
where the shape of minimum energy is unique and is a prolate 
as pointed out in Sec. IIB. 

We compare the results for the uniformly-charged shells 
with conducting spheroidal shells. The energy expression for 
the volume-constrained conducting shell problem is the same 
as Eq. (25) except that the functions a and c are now deter- 
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mined by Eqs. (30) and (27) respectively. Similar to Fig. 2, 
we find the Coulomb energy of a conducting shell is always 
lower than that of the homogeneously charged shell. Again, 
we attribute this to the fact that for a conducting shell, the 
charges will move until the shell surface becomes an equipo- 
tential, and this movement always lowers the Coulomb en¬ 
ergy. Finally, we note that in Fig. 3, the energy profiles for 
the two different scenarios (homogeneous and conducting) are 
similar in shape which is in contrast with Fig. 2 for the area- 
constrained case. 


III. ENERGY ANALYSIS OF PERTURBED SPHERES 


Results of the preceding sections indicate that a uniformly- 
charged sphere is locally stable to perturbations towards a pro¬ 
late or oblate spheroid if the deformations preserve the surface 
area. The stability of the charged sphere for arbitrary pertur¬ 
bations in shape that preserve the area and maintain the unifor¬ 
mity of the surface charge distribution is not addressed yet. In 
this section, we explore this problem by deriving the expres¬ 
sion for the energy variation of the shell due to small generic 
perturbations on the spherical shape. 

Consider a sphere of radius Rq that is represented by 

Ro{9,(j)) = {Rosin9cos<j),Rosin9sincj),RoCOs9), (33) 

where 9 G [0, tt] and (p G [0, 27r). A generic small perturba¬ 
tion h{9, (p) on the sphere produces a shape that can be repre¬ 
sented by 

R{9, P) = Ro{9, p) + h{9, p)N{9, p), (34) 

where N{9, p) = Rq^Rq{ 9, p) is the unit normal vector on 
the sphere. Note that R{9,p) = Rq{ 9,P){1 + h{9,p)/RQ) 
and we consider perturbations to be small when they satisfy 
the condition |/i|/i?o ^ 1- The metric tensor of the unde¬ 
formed sphere is 


with its determinant go = ( 50 ) 11 ( 50)22 = Rq sin^ 9. The de¬ 
terminant of the metric tensor of the deformed shape defined 
inFq.(34) is [14]: 


9 = 90 + Sg, 


(36) 


where 

6g = go [4hH + {goY^h,h, + + 2K)\ + 0{h^), 

(37) 

with H = l/i?o and K = (l/i?o)^- 

We assume that the charges on the surface interact via a 
pair potential V{r) which only depends on the distance r 
between them. We consider a homogeneous surface and al¬ 
low only those deformations that keep the area unchanged. 
Thus, the charge density is the same even in the deformed 
shape. The charge elements on the unperturbed shape are 


dqo = (JodAo = ao^/god^x where ao = Q/Ao is the uni¬ 
form charge density. The charge elements on the perturbed 
shape are dq = aodA = uoy/gd^x. Note that d^x is a short¬ 
hand for d9dp. The change in the interaction energy brought 
about by the deformation of the charged sphere is 


2 n 

H[h] = yJ \/9{x)d?xyJg{x')d^x' V R(x') - R(x) 

2 n 

- ^ / \/go{x)d^x\Jgo(,x')d^X V Ro(x) - Ro(x) 


(38) 


where the integration is over xi = 9, X 2 = p, x\ = 9', x '2 = 
p'. We note that if 77 <0, the energy of the perturbed sphere 
is lower than that of the spherical system. The distance be¬ 
tween two arbitrary points on the deformed sphere is 


R{x') — R{x) 

(^Ro{x') — Ro{x)^ + (ji{x')N{x') — h{x)N ix)^ 


= |ro -I- 5r| = ro -I- Ar, 


(39) 


where tq = Ro{x') — Ro{x) and 5r = h{x')N{x') — 
h{x)N{x). The r in Fq. (39) can be expanded in terms of 
5r according to 


f + 5f 


/ + A/, 


(40) 


where / = | / | and 


A/ 


/■ 6f 1 Sf- Sf 
f 2 f 


1 if- 0? 

2 /3 


+ 0{5f). 


(41) 


Therefore, 


Ar = {h{x) + h{x')) + (h{x) - h{x')f (42) 

2 7?o 2ro 

+ ^ (h(x) + h(x')f + O(h^). 

Zd -iLQ OjCq 

We expand V (r) and ^/g{¥)g{P^ in terms of h up to the 
quadratic order: 


V(r) = V (ro + Ar) 


(43) 

C>(Ar3), 


'/^Vdix') = Vg^o X (44) 

2 Uo 5o 9o9o) sUo 9o) 

where the prime on functions go and 5g denotes that the vari¬ 
able of the function is x', e.g., g'o is short for go{x'). Note that 

6g contains terms linear in h. _ _ 

We formally write \/ 9 {x)y/g{x') = 

\/ 3 o 5 o(l + a(/i)-b/3(/i2)) and U(r) = V{ro) + A{h) + 
B{h?), where a{h) and A(h) are terms linear in h and /3(h^) 


4-0(h3), 
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and B{h?) represent terms quadratic in h. Eq. (38) can 
therefore be written as 

2 n 

H = Y J cPxcPx'^/gogo [A{h) +a{h)V{ro)] 

+ Y / d^xif^'Jgog'o + a{h)A{h) + P{h'^)V (ro)] 

+Oih^). (45) 


Using Eqs. (37) and (42), we obtain A{h) = 2 ^ + 

h') and a{h) = -^{h + h'). Up to the linear term in h, 
Eq. (45) becomes 


2 


d9d(j)d9' d4> sin 9 sin 9’ 


2U(ro) + 


rp dV{ro) 
2 dxQ 


x{h{9,cj)) + h{9\^'))+0{h^). (46) 


We find that app is always negative, implying a uniform 
shrinking of the shell to preserve the area. When only uni¬ 
form shrinking (or expansion) is allowed, aio = 0 for I > 
0. Eor this case, from Eq. (50), we obtain the relation: 
2-\/4^appi?p + Upg = 0. This equation has two solutions, 
opp = —A^/ttRq, which is unphysical as |app| > i?g, and 
apo = 0. The latter solution indicates no deformations on the 
sphere. 

Now we consider perturbations characterized by two modes 
Opp and flip: 


h{9,(l)) =+ aio\ cos 9, (51) 

i/47r V 47r 

where 9 S [0,7r]. We note that the modes ago and oip are 
coupled due to the constraint of fixed area, Eq. (50), leading 
to the relation: 


Note: 

ro(9,^,9',<^')= (47) 

BovW — cos 9 cos 9' — cos(()) — Ip') sin 9 sin 9'\. 

The existence of terms linear in h in Eq. (46) suggests that 
the spherical shape is not necessarily an energy extreme, open¬ 
ing the possibility of shapes with lower energy than that of the 
sphere. Eor V{rg) = l/rg, the term in the square bracket 
becomes 


-Ogp + 2\/TrRoaoo + afg — 0. (52) 

Recall that agg < 0 and the above equation implies aip can 
assume positive or negative values. Physically, the perturba¬ 
tion h of Eq. (51) corresponds to a buckling of the spherical 
shape. If we identify 6* = 0 as the north pole and consequently 
0 = TT as the south pole, we find that the north pole is buckled 
inward for aip < 0 and the south pole is buckled inwards for 

oip > 0. 


2U(rp) -b irg 


dVjrg) 

drg 



(48) 


We find that for a = 1, which represents the Coulomb poten¬ 
tial, this term is non-zero. This indicates that the spherical sys¬ 
tem may be unstable to the long-range Coulomb interaction. 
We note that for a = 4, the above term vanishes implying 
there are no linear terms in h for this potential. Furthermore, 
the sign of the energy changes as the value of a goes beyond 4, 
implying the distinct stability behaviors of the spherical shape 
for long- and short-range potentials. 

We now investigate if there exists a perturbation that can 
lower the Coulomb energy of the uniformly-charged sphere, 
in other words, a deformation for which H < 0. We consider 
a generic form for the perturbation h represented hy h{9,p) = 
X]; m where are the spherical harmonic 

functions with Z = 0,1, 2,... and m = —I, —l+l,... ,1. Here 
aim are the unknown expansion coefficients or modes. As h 
is taken to be a small perturbation, we require \aim\ ^ Rq- 
We take m = 0, thus examining axisymmetric deformations. 
The constraint of fixed area leads to the relation between the 
coefficients aig. The variation of area [14] is 


SA — 2^/^aggRg + 

Im ^ ' 

Setting 8A = 0 leads to 

2^/^aggRg = — lUimP f ^ 4” 2^^^ 

Im ^ ^ 



FIG. 4. The total Coulomb energy change H vs the coefficient ago- 
The inset shows the variation of R with aip. The shape of the per¬ 
turbed sphere is characterized by the two modes ago and aio, which 
are coupled due to the constraint of the fixed area. 


The change in the total Coulomb energy H of the system 
(up to linear terms in h) for this perturbation can be calculated 
from Eq. (46) and we find: 


_ 3^0 i^o 
2^4^ 


d9dpd9' dp' sin 9 sin 9 


, opp -F oiov/ScosS 




(53) 

We note that the shape of the perturbed sphere is independent 
of the sign of the coefficient aip; to reverse the sign of aip 
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is to rotate the shape by tt. This suggests that terms linear in 
Cio should be absent in H. By making use of the parity of the 
integrand in the four quadrants: {0 G [O,7r/2],0' £ [0,7r/2]}, 
{9 G [0,7r/2],6»' G [7r/2,7r]}, {9 G [7r/2, tt], 6>' G [0,7r/2]}, 
and {9 G [7r/2,7r],0' G [7r/2,7r]}, we find that the integral 
involving the oio cos 9 term in Eq. (53) vanishes. Noting that 
the remaining integral (involving the ago term) in Eq. (53) can 
be read as the total Coulomb energy of a uniformly-charged 
spherical shell, we arrive at the following analytical result for 
H: 


H = (54) 


and note that and dH have the same sign and differ by 
a constant prefactor. Upon substituting h from Eq. (51) in 
Eq. (57), the latter becomes 




ooo 


73 


Olo 


(58) 


We follow a similar procedure to obtain the change in the lo¬ 
cal Coulomb energy density Hs at the south pole. We find 
that Hs is given by Eq. (58) with aio replaced by —uiq. Ei- 
nally, noting that the two modes ago and oio are coupled via 
Eq. (50), we express Hjy and Hs as functions of aoo: 


Because ago < 0, we see from Eq. (54) that iJ < 0. In other 
words, the perturbation in Eq. (51) lowers the Coulomb en¬ 
ergy of the original unperturbed spherical shape. In Eig. 4, we 
plot H vs agg for cTg = 1 and i?o = 1. Recalling that agg 
characterizes the amount of buckling in the deformed shape, 
we find that the system energy keeps decreasing with the rise 
in the buckling of the shape. Noting that agg and aw are cou¬ 
pled via Eq. (50), we obtain H as a function of aio using 
Eq. (54): 

H = 127/ (^^4TTRl-2al^ - sj^ ttRI^ ■ (55) 

As expected, H is found to be independent of the sign of uiq. 
In the inset of Eig. 4, we show the variation of H with uiq. 
The above calculations demonstrate that a uniformly-charged 
sphere is electrostatically unstable to a buckling-type defor¬ 
mation in the constraint of fixed area. 

Alongside the computation of H, which quantifies the 
change in the total Coulomb energy, it is instructive to exam¬ 
ine how the Coulomb energy changes locally at specific points 
on the shell surface as a result of the deformation proposed in 
Eq. (51). We define the local Coulomb energy at a point as the 
interaction energy of the charge element at that point with all 
the other charges on the surface. Eor the uniformly-charged 
sphere (undeformed state), the local Coulomb energy is the 
same at all points on the surface. 

We discuss, without any loss of generality, the case of 
oio < 0 for which the north pole {9 = 0) is buckled inwards. 
Eor simplicity, we perform calculations for two points on the 
surface: the north and the south poles. Eor the perturbation h 
given by Eq. (51), we find the change in the local Coulomb 
energy (relative to the original spherical conformation) of a 
charge element located at the north pole to be: 

dH = -(jgdA' [ Rg sin 9d9d(j)-^ {h {9, (j)) + h (0, 0)), 

2 J 2rg 

(56) 

where dA' = 7 9oixN)d‘^XN is the area element associated 
with the charge in the undeformed (spherical) conformation 
and rg = 2Rgsm{9/2). We are primarily interested in the 
sign of dH', if this interaction energy is negative (positive), 
that would imply the local Coulomb energy for the deformed 
shape is lowered (raised) relative to the original sphere. We 
define 

Hn = cfgRg I d9d(j)sm9-^{h{9,(j)) + h{0,0)) (57) 

J 2rg 


Hn,s = Gy/nal ^aoo T (-47^aooi?o - floo)^ > 

(59) 

where the — (-f) sign corresponds to the north (south) pole. 
We plot iJjv and Hs vs agg in Eig. 5 for o-q = l,Rg = 1. 
We find that as the magnitude agg of the deformation (buck¬ 
ling) is increased, the local Coulomb energy is lowered at the 
north pole (Hjv < 0). On the other hand, at the south pole, 
the local Coulomb energy is higher relative to its value in the 
undeformed spherical case (Hs > 0). 



®oo 

FIG. 5. Hn (dashed green) and Hs (solid red) vs aoo for the pertur¬ 
bation h characterized by two modes aoo and aio. The above plot is 
for the case when aio < 0, which corresponds to the inward buckling 
of the north pole. See text for the meaning of the symbols. 


IV. CHARGE RENORMALIZATION IN SPHEROIDAL 
SHELLS 

In previous sections, we analyzed the Coulomb energy 
of uniformly-charged spheroidal and sphere-like shells and 
determined the conformations that correspond to the lowest 
Coulomb energy under a given geometric constraint. The 
shell was considered to be an isolated system in vacuum with 
charges embedded on its surface. In realistic settings, how¬ 
ever, we expect the shells to be in an environment where 
the surrounding medium (solvent) contains counterions that 
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neutralize the shell charge, rendering the overall system elec¬ 
troneutral. In this light, results obtained thus far are reliable 
in the event that the associated counterions remain in the bulk, 
far from the shell surface. In situations where a significant 
fraction of the total number of counterions condense on the 
shell, the free energy of the shell-counterion system, deter¬ 
mines the equilibrium shell conformations. We note that at 
infinite dilution, in the spherical case, the entropy is expected 
to dominate the shell-counterion Coulomb attraction, leading 
to no counterion condensation on the shell surface. However, 
at finite shell concentrations (volume fractions), condensa¬ 
tion is expected to occur, even in salt-free settings [15-17]. 
The condensation of the counterions can be viewed as renor¬ 
malizing the (bare) charge on the shell [15]. This renormal¬ 
ized charge and consequently the behavior of the equilibrium 
free energy of the shell-counterion system can be obtained 
in a qualitative way by using the Manning-Oosawa two-state 
model [ 8 , 9, 16, 18]. In this section, we use this two-state 
model approximation and compute the renormalized charge 
on uniformly-charged oblate and prolate spheroidal shells at 
finite shell concentrations in salt-free settings and find the 
variation of the equilibrium free energy of the system as a 
function of the shell aspect ratio. 

We investigate the effects of counterion condensation on 
spheroidal shells of equal area and spheroidal shells of equal 
volume. For either case, we take the area (volume) to be 
constrained to that of a sphere of radius R. We consider a 
Wigner-Seitz (WS) cell of volume Vws containing a single 
shell of volume H, with Q charge on its surface, placed at 
the center. We work with finite values of the shell concentra¬ 
tion 77 = fl/Vws- The cell also contains N counterions, each 
of charge Q/N making the overall shell-counterion system 
electroneutral. The counterions are separated into two dis¬ 
tinct groups: free ions and condensed ions. The condensed 
counterions are restricted to have translational motion in a 
thin layer of volume 14 = ^(A, R)h surrounding the shell, 
where A(A, R) is the area of the shell and b is the thickness of 
the layer. Note that as we restrict our analysis to shells con¬ 
strained to a fixed area or volume, the area A (and the volume 
H) of the shell can be considered as a function of the aspect 
ratio A and R. Free ions occupy the available space in the 
WS cell which in the dilute limit can be approximated to be 
the volume of the cell. We choose experimentally relevant pa¬ 
rameters: total charge Q = 600 electron units (which amounts 
to ~ 100 mV of surface potential), N = 1000 counterions, and 
i? = 10 nm. Calculations are performed for room tempera¬ 
ture T = 300 K and we take water as the dielectric medium 
surrounding the shell. 

Let a be the fraction of counterions that condense. Clearly, 
{1 —a)N ions remain free in the bulk. Further, the condensed 
ions neutralize the surface charge on the shell reducing the 
net charge to (1 — a)Q. We approximate the WS cell to be 
spherical with volume Vws = ( 4 / 3 ) 7 ri?^s, with i?ws being 
the radius of the cell that gets determined by the shell volume 
fraction (concentration) rj. We vary the shell concentration rj 
from 10“^^ to 10“^. We write the free energy (in units of 


ksT) associated with the shell as: 


F{a, A) = (1 - afU{X) + aN\n 
-f (1 - a)N In 


aNA^ 


A{\R)b 
{l-a)NA^'\^ 


- aN 




WS 


-il-a)N, 


(60) 


where A is the thermal de Broglie wavelength. The first term 
is the electrostatic potential energy of the shell with renormal¬ 
ized charge (1 — a)Q. The function U (shown below) rep¬ 
resents the bare (unrenormalized) Coulomb energy and is de¬ 
termined based on particular geometric constraint employed. 
The second and third terms stem from the entropic contribu¬ 
tion of the aN condensed ions, and the last two terms cor¬ 
respond to the entropy of (1 — a)N free counterions. Note 
that within this model, the entropy of both free and condensed 
ions is assumed to be that of an ideal gas. Also, just like U, 
the form of the area A of the shell and its volume fl depends 
on the constraint applied as shown below. 

For the case of shells subject to the equal area constraint, 
the function 17(A) reads 


C/(A) 



1b17o,a('\/1 — A^, R,Q) 0 < A < 1 
lBC/p,A(yA^^/A,i?,Q) A>1, 

( 61 ) 


where Uo^a and Up^A are available from Eq. (24), which pro¬ 
vides the expression for the Coulomb energy of uniformly- 
charged spheroidal shells of equal area, and the above re¬ 
sult is expressed in units of kgT by introducing the Bjerrum 
length /fl. Recall that the values of A > 1 correspond to pro¬ 
late spheroids and 0 < A < 1 region corresponds to oblate 
spheroids. Owing to the constraint, the area of the shell is 
simply A(A, R) = AnR?. The shell volume H follows from 
Eqs. (26), (19), (29), and (22): 


H(A,i?) 


Oo(yrT7^,a(yrTrA2,i?)) 0<A<1 
Op (VA2 - 1/A, c (VA2 - 1/A, R)) A>1. 

(62) 


Eor the volume-constrained problem, U (A) follows from 
Eq. (32): 


t/(A) 



IbUo,v{VT^, R,Q) 0 < a < 1 

IbUp,v(.V>^^/\,R:Q) A>1. 


(63) 

Eollowing Eqs. (7), (12), (27), and (30), the area function A 
for this case becomes 


( <{VT^,a{VT^,R)) 0<A<1 

^ ’ '' t (VA2 - 1/A, c (yFTTi/A, R)) A>1. 

(64) 

Einally, owing to the constraint, the shell volume is simply 
n = AttR^/S. 

We approximate the thickness b of the condensed layer by 
the Gouy-Chapman (GC) length b = 1/{2'!tIb(j), where a is 
the unrenormalized charge density on the shell surface [18]. 
Higher charge density or longer Bjerrum length leads to a 
stronger shell-counterion attraction implying a thin condensed 
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layer; this is indeed reflected when b is chosen as the layer 
thickness as seen from the above expression. We also note 
that the GC length is a length scale associated with the planar 
interface and hence our analysis is limited to the regime where 
b is shorter than the characteristic lengths associated with the 
shell. We have carried out the following analysis by choosing 
the Bjerrum length /b as the thickness of the condensed-layer 
and we And no changes in the conclusions reached below. 

The free energy F in Eq. (60) can be considered as a func¬ 
tion of A and a. For a given A (shape), we minimize the free 
energy with respect to a to And the fraction of counterions that 
condense on the shell. We obtain the extremum condition: 


—.^(1 — a) -I- Win 


/ g 1 Vl{X,R) \ 

\ 1 — a ry A(A, R)b J 


= 0 , 


(65) 


where ^ = 2U{\) measures the strength of the Coulomb in¬ 
teractions and p is the volume fraction of the shells given by 



( 66 ) 


For a given A and rj, we solve Eq. (65) using Mathematica 
and obtain a as a function of A. We carry out the study for 
a wide range of shell volume fractions ranging from 10 “^^ 
to 10“"^. Using the value of a, the renormalized electrostatic 
energy U of the shell at equilibrium is known from Eq. (61) or 
Eq. (63) (depending on the constrained problem under study), 
by replacing Q with (1 — a)Q'. 

U{\Q) = U{\{l-a)Q). (67) 


Employing the above result and the equilibrium value of the 
condensate fraction a (obtained as the solution of Eq. (65)), it 
is easy to show that the difference in the equilibrium free en¬ 
ergies of a spheroidal shell and a spherical shell, dF, defined 
as 


dF{X) = F{X) - F(A 1) ( 68 ) 


is given by: 




I — a 


1 — a. 


U, + Wln-^ “ 


1 — as 


(69) 


In Eq. (69), as and Us denote, respectively, the values of the 
condensate fraction and the renormalized Coulomb energy for 
a spherical shell. 

We now analyze the variation of dF as A is changed. We 
consider A values from 0 to 4 like in the study of the Coulomb 
energy of isolated, homogeneously-charged spheroidal shells 
recorded in Figs. 2 and 3. We begin with the case of fixed 
shell area. In Fig. 6 , we plot dF, computed from Eq. (69), as 
a function of A for various values of rj. The red solid line is 
the Coulomb energy dU of isolated spheroidal shells of equal 
area measured relative to the Coulomb energy of the spherical 
shell with identical parameters. This line acts as a reference 
curve to all other lines which are the result of taking ion con¬ 
densation into consideration. Recall that for the counterion- 
free case, the spherical shape is a local energy minimum as 



FIG. 6 . Equilibrium free energy difference, dF, between a 
homogeneously-charged spheroidal shell and a sphere of identical 
parameters (z = 0.6, N = 1000, 7? = 10 nm) as a function of the 
aspect ratio A for the area-constrained system. The effects of charge 
renormalization due to counterion condensation are included via a 
two-state model analysis calculation. The red curve is the case for 
the counterion-free system and acts as a reference line. All other 
curves take into account ion condensation on the shell surface and 
correspond to different values of the shell volume fraction p. Results 
are shown for 77 = 10 ~^^ (green squares), 10 “^'’ (blue circles), 10 “® 
(orange triangles), 10 “® (cyan inverted triangles), and lO”"^ (brown 
diamonds). We find that the spherical conformation, which is a lo¬ 
cal minimum for the isolated spheroidal system (red line), becomes 
a free-energy maximum at 77 = 10“^. 


evidenced by the dU plot. For low rj (= 10“^^), which cor¬ 
responds to a very dilute system, we find that the dF curve 
(green squares) lies in the vicinity of the no-condensation re¬ 
sult {dU). However, as 77 rises, we observe significant devi¬ 
ations from the unrenormalized energy curve for both oblate 
(A < 1) and prolate A > 1 regions. 

We find that for all shapes (A), as the volume fraction 77 
is increased, a increases, that is more counterions condense 
on the shell surface. A major consequence of the enhanced 
charge renormalization is the reduction of the positive free en¬ 
ergy difference between the spheroidal shell and the sphere 
(see Fig. 6 inset). Further, for high 77 values (77 = 10“®, 
10 “"^), we find that all prolate shells have less free energy as 
compared to the sphere which is in stark contrast from the 
no-condensation result (red line). For the same 77 values, we 
also find that oblate shells with small aspect ratios have lower 
free energy than a spherical shell {dF < 0). For the volume 
fraction of 77 = 10 “'^ (brown diamonds), the spherical shell 
has the maximum equilibrium free energy among all shapes. 
Thus, according to the above analysis based on the two-state 
model, for the area-constrained shell system, counterion con¬ 
densation has a profound effect in modifying the energy land¬ 
scape associated with the isolated shell, favoring the forma¬ 
tion of spheroidal structures over the spherically-shaped ones 
as the shell volume fraction is increased. 

Fig. 7 shows the variation of dF with A for different 77 val¬ 
ues in the case of fixed-volume constraint. Once again, the 
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FIG. 7. Equilibrium free energy difference, dF, between a 
homogeneously-charged spheroidal shell and a sphere of identical 
parameters (z = 0.6, N = 1000, 7? = 10 nm) as a function of 
the aspect ratio A for the volume-constrained system. The effects of 
charge renormalization due to counterion condensation are included 
via a two-state model analysis calculation. The red curve is the case 
for the counterion-free system and acts as a reference line. All other 
curves take into account ion condensation on the shell surface and 
correspond to different values of the shell volume fraction rj. Results 
are shown for ?7 = 10 “^^ (green squares), 10 “^° (blue circles), 10 “® 
(orange triangles), 10 ~® (cyan inverted triangles), and lO”"^ (brown 
diamonds). We hnd that in the event of ion condensation, for all val¬ 
ues of 77 , the spherical conformation continues to have the highest 
free energy among all spheroidal shapes. 


solid red line corresponds to the counterion-free system and is 
the Coulomb energy of the uniformly-charged spheroidal shell 
measured relative to the electrostatic energy of the sphere. 
As Fig. 3 shows, the spherical shape is the conformation of 
maximum energy for this particular constraint. This conclu¬ 
sion remains unchanged when we include the effects of charge 
renormalization via the two-state model analysis as seen from 
Fig. 7. For a very dilute system, -q = 10“^^ (green squares), 
we find that the dF curve lies in the vicinity of the dU line. 
Increasing 77 leads to a rise in a which is seen to weaken the 
(negative) difference between the equilibrium free energies of 
spheroidal and spherical shells (see Fig. 7 inset). Thus, judg¬ 
ing by the variation of dF determined by the two-state model 
analysis, we find that oblate and prolate shaped structures con¬ 
tinue to be energetically favored over the spherical confor¬ 
mation in the event of counterion condensation. Further, the 
monotonic trend of free energy decrease with increasing the 
eccentricity of the spheroidal shell observed for the isolated 
shell system is seen to persist in the wake of counterion con¬ 
densation as well. 

We note that as the counterions are mobile, under certain 
conditions, the charged shell is better approximated as a con¬ 
ducting surface as opposed to a homogeneously-charged one. 
In that event, the above analysis can be carried out using the 
expressions of electrostatic energy for conducting spheroidal 
shells provided in Eq. (25). We observe that for all values 
of 77 , the conducting spheroidal shell has a lower free energy 


upon deformation in comparison with the homogeneously- 
charged shell (see Ref. 7 for details of this calculation for 
oblate spheroidal shells under the constraint of constant vol¬ 
ume). Further, we find that the main conclusions regarding the 
effects of counterion condensation reached above for either 
constraints remain unchanged when we repeat the two-state 
model analysis assuming that the shell is an equipotential sur¬ 
face. 

Finally, we note that the study of condensation effects 
based on the above two-state model employs a number of ap¬ 
proximations. For example, the free energy associated with 
this model, Eq. (60), does not take into account the shell- 
counterion and counterion-counterion Coulomb interactions 
explicitly. Another approximation is to put the shell in a 
counterion-only, salt-free environment. Also, the distribution 
of counterions around the spheroidal shell is considered to be 
isotropic which is clearly a simplification for shapes that devi¬ 
ate significantly from the spherical conformation. In this light, 
we view the above results as qualitative. Quantitative results 
that address many of the aforesaid simplifications can be ob¬ 
tained by employing approaches based on the solution of the 
Poisson-Boltzmann equation for spheroidal geometry [19]. 


V. CONCLUSION 

We report the exact expression for the electrostatic energy 
of a uniformly-charged spheroidal shell. We analyze the vari¬ 
ation in the electrostatic energy as the aspect ratio of the 
shell is changed from 0 to 00 for the area-constrained and 
volume-constrained cases. The prolate spheroidal shell with 
its major-axis length stretched to infinity is found to have 
the lowest Coulomb energy among spheroidal shells of equal 
area. Further, we reveal the non-monotonous variation in the 
Coulomb energy when a spherical shell is elongated to a pro¬ 
late spheroid keeping the shell area fixed. For spheroidal 
shells that have the same volume, a sphere has the highest 
Coulomb energy. In addition, our perturbation calculations 
show that there exist area-conserving buckling-type deforma¬ 
tions on the sphere that can lower the total Coulomb energy. 
For the spheroidal shell system, we use a two-state model of 
free and condensed ions to evaluate the renormalization of the 
shell charge due to counterion condensation. We find that ion 
condensation has a significant effect in modifying the free en¬ 
ergy landscape with spheroidal structures being favored over 
a sphere of the same area as the shell volume fraction is in¬ 
creased. These results add to the theoretical foundation re¬ 
quired to understand the control of spheroidal shapes in mate¬ 
rials using electrostatics in combination with other forces such 
as those arising due to the elastic nature of the material [7]. 
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Appendix A: Coulomb energy of uniformly-charged prolate 
spheroidal shells 

In this section we derive the expression shown in Eq. (9) 
for the Coulomb energy of a homogeneously-charged prolate 
spheroidal shell. We consider prolate spheroidal coordinates 
M, V, (j), which are related to the Cartesian coordinates x, y, z 
by 

a; = cCp sinh(u)sin(?;)cos((/)), (Al) 

y = cCp sinh(u)sin(?;)sin((/)), (A2) 

2 ; = cepCosh(u)cos(?;), (A3) 

where 

0 < U < 00 , 0 < V < TT, —TT < (j) < TT. (A4) 

The set {u^v,(j)) uniquely characterizes a point in the 3- 
dimensional space. It is straightforward to show that the met¬ 
ric coefficients are 

hu = hy = cCp \/ sinh^M + sin^u, hcj, = cCp sinhu sinu 

(A5) 

using which the form for the Laplacian is readily ob¬ 
tained to be [ 12 ] 


A closer examination of these equations reveals the general 
solution for the potential to be; 

00 

^{u,v) = y^{AnPn{co&hu) -f (coshu))P„(cosu) 

n—0 

(AlO) 

where n is an integer, and Qn are Legendre functions of 
the first and second kind respectively, and A„ and are un¬ 
known coefficients. In order to ensure that the solutions are 
bounded in the interior and exterior regions of the spheroid, 
we find that A„ must vanish in the domain u > uq and 
Bn = 0 in the region where 0 < u < uq. We thus have the 
following form for the potential inside and outside the oblate 
shell; 

00 

T>in(u,'!;) = A„P„(coshu)P„(cosu), (All) 

n—0 

00 

$out(u, u) = y] P„Q„(coshit)P„(cosu). (A12) 

n—0 

The potential must be continuous at the shell surface u — 
Uq = 0, that is, $in(uo,u) = 4’out(uo,u). This boundary 
condition leads to the relation 

A„P„(coshito) = P„Q„(coshuo) (A13) 

forn = 0,1, 2,.... Note that coshuo = l/cp. The discontinu¬ 
ity in the normal component of the electric field at the charged 
surface provides another boundary condition; 


= 


c 2 ep(sinh^it -|- sin^u) 


Ml 

\ sinhu Su 




1 < 52 $ 

c^CpSinh^u sin^u 


1 S 
sinu Sv 



(A 6 ) 


The prolate spheroidal shell in these coordinates is repre¬ 
sented by the simple equation u = uq, where uq is connected 
to the eccentricity via the relation 


sechuo = Cp. (A7) 

The region of space interior to the spheroid corresponds to 
the values 0 < u < uq and the exterior region is repre¬ 
sented by the u > uq domain. We begin by finding the elec¬ 
trostatic potential generated by the uniformly-charged prolate 
shell represented by the equation u = uq. Since there is ax¬ 
ial symmetry in the problem, the electrostatic potential will 
depend only on coordinates u and v. Writing the solution 
as = U{u)V{v) and substituting it in the Laplace 

equation = 0, we find, upon using Eq. (A 6 ), that the 
variables separate and the functions U and V satisfy the dif¬ 
ferential equations; 

1 s / su\ 

—— — sinhu— -n[/ = 0, (A 8 ) 

sinhit ou \ ou J 

1 S ( (514\ 

^—-— I sinu— I -I- nV = 0. (A9) 

sinu dv \ ov ) 


— U • V$out + U • V$in = dTTCr, atM = Mo- (A14) 


Using the expression for the gradient in prolate spheroidal co¬ 
ordinates, the above equation becomes 


1 

hu Su 


U — Uq 


1 S^o 


hu Su 


= 47ra. 


(A15) 


Employing equations (All), (A12), and (A5), the above 
boundary condition leads to the relation 


^(AM(l/ep) 

n—0 


BnQ'nO/ep)) Pn(cOSv) 


47r creep 



1 — BpCOS^U. 


(Alb) 


Eqs. (A13) and (Alb) allow us to evaluate the undetermined 
functions A„ and Bn- We first eliminate in favor of A„ 
using Eq. (A13) obtaining 


AnPn(l/ep) 

f^n(l/ep) 


(A17) 


Substituting Bn from above in Eq. (Alb), we obtain 


n—0 


AM(l/ep) 


AnPnO/ep) 

f^n(l/ep) 


QnO/ep) 


Pn(cOSv) 


dTTcraep 



1 — e^cos^ 


(A18) 
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Using the fact that the Wronskian of (z) and Qn (^) is given Appendix B: Coulomb energy of uniformly-charged oblate 
by spheroidal shells 


lU(P„(z),Q„(z)) = -^, (A19) 

1 — 

Eq. (A18) can be simplified to 


AnPnjcOSv) 

n=0 


dTTCTC 



Cp 



COS^ti. 


(A20) 


Multiplying both sides of the above equation with 
Pi (cosu)sinu and using the orthogonality relation; 

r 2 

/ Pn{cosv)Pi{cosv)sinvdv = - — —Snl, (A21) 

Jo 2n+l 


we obtain 


. 2n+l 
Ari = —-— Anac- 


Q„(l/ep)iJ„(ep), (A22) 


where Hn{ep) is the integral 

Hn{ep) = J Y^l — CpCos^u P„(cosz;)sinu dv. (A23) 


It is easily checked that for odd n the integral in the above 
equation vanishes. Hence, we have Ai = ... = Q. 

Using An in Eq. (A17), Bn is known as well, and conse¬ 
quently from Eqs. (All) and (A12), we obtain the desired 
electrostatic potential at any point in space. 

Eor the computation of the electrostatic energy, the knowl¬ 
edge of the potential on the shell surface suffices. Using 
Eqs. (All) and (A22), we obtain the surface potential as 


dTrcrcW 1 — Cp 

^sheii(w, ep,c) =-X 

ep 

^n(l/ep)Qn(l/ep)i?„(ep)P„(cosu), 

n^even 

(A24) 


where the summation is over even integers n = 0, 2 ,4,.... 
The electrostatic energy of a charged spheroidal shell can be 
written as 


C/= (A25) 

The shell surface area element in the prolate spheroidal coor¬ 
dinates is given by 


dA = hyh^dvdj) = — — CpCos^u sinu dvdj) 

(A26) 

where the second equality follows from Eq. (A5). After sub¬ 
stituting this expression for the area element in Eq. (A25), us¬ 
ing Eq. (A24), and changing the variable from a to Q, we 
obtain the expression in Eq. (9). 


We now derive the expression for the electrostatic energy of 
a uniformly-charged oblate spheroidal shell. As most of the 
steps involved in this derivation are analogous to the above 
derivation for the prolate case, we will present only the key 
steps of the procedure. We begin by employing the oblate 
spheroidal coordinates u, v, 4>, which are related to the Carte¬ 
sian coordinates x,y,zhy 

a; = aCo cosh(M)sin(v)cos(())), (Bl) 

?/= aCo cosh(M)sin(v)sin((^), (B2) 

z = aCo sinh(M)cos(v), (B3) 

where 

0 < tt < OO, 0 < V < TT, —TT < (j) < TT. (B4) 

The set {u^v,(j)) uniquely characterizes a point in the 3- 
dimensional space. It is straightforward to show that the met¬ 
ric coefficients are 

hu = hy = aCo V sinh^M -f cos^u, = aCo coshu sinu 

(B5) 

using which the form for the Laplacian = 0 is readily 
obtained. The oblate spheroidal shell in these coordinates is 
represented by the simple equation u = uq, where uq is con¬ 
nected to the eccentricity via the relation 

sechuQ = Cq. (B6) 

The region of space interior to the spheroid corresponds to the 
values 0 < u < itQ and the exterior region is represented by 
the u > Uq domain. 

Once again, we start by determining the electrostatic poten¬ 
tial generated by this uniformly-charged oblate shell. Since 
this system has axial symmetry, the electrostatic potential cre¬ 
ated by the oblate spheroid will depend only on the coordi¬ 
nates u and V. Employing separation of variables we can write 
the solution as <I>(m, u) = C/(u)U(t;), upon which the Laplace 
equation separates into two differential equations for u and v. 
A closer examination of these equations reveals the general 
solution for the potential to be: 

C30 

4>(u,u) = y^(A„P„(isinhu) -f sinhu))P„(cosu), 

n—0 

(B7) 

where Pn and Qn are Legendre functions of the first and sec¬ 
ond kind respectively, and An and Bn are unknown coeffi¬ 
cients. In order to ensure that the solutions are bounded in 
the interior and exterior regions of the spheroid, we find that 
An = 0 in the domain u > uq, and Bn = 0 when 0 < rt < uq. 
We thus have the following form for the potential inside and 
outside the oblate shell: 

OO 

$in = y] AnPn{i sinhu)P„(cosu), (B 8 ) 

n—0 

OO 

^OUt - BnQn{i sinhu)Pji(cosu). (B9) 

n=0 
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The boundary condition that the potential must be continu¬ 
ous at the shell surface u — Uq = 0 leads to the relation 

AnPn{i sinhito) = BnQn{i sinhito) (BIO) 


Using Eqs. (B14) and (B15) in Eq. (A25), and changing the 
variable from a to Q, we obtain the energy expression in 
Eq. (14). 


for n = 0,1,2,.... Note that sinhrto = (l/eo)-\/l — e^. 
The discontinuity in the normal component of the gradient of 
the electric potential at the charged surface provides another 
boundary condition which upon employing the expression for 
the gradient in oblate coordinates becomes: 


OO 

{AnPn(i sinhwo) — BnQni'^ sinhuo)) icoshuoPTi(cosu) 

n—0 


= Airaa 



eoSin^u. 


(BID 


Using Eqs. (BIO) in (Bll), we can eliminate Bn in favor of 
An and solve for the latter, obtaining 


A 


n 


2n + 1 Airaa i 

2 Co 


Qn 



Dl (t^o) ; 


(B12) 


where In{eo) is the integral 


In (t'o) 




eoSin^?;P„(cosu)sinz; dv. 


(B13) 


As before, in arriving at this result we employed the prop¬ 
erty of the Wronskian of the Legendre polynomials and their 
orthogonality relation. It is easily checked that for odd n, 
In{eo) vanishes, implying Ai = A 3 = A 5 ... = 0. Using A„ 
in Eq. (BIO), can be evaluated as well and consequently 
from Eqs. (B 8 ) and (B9), the electrostatic potential is known 
everywhere in space. 

Eor the computation of the electrostatic energy, the knowl¬ 
edge of the potential on the shell surface suffices. Using 
Eqs. (B 8 ) and (B12), we obtain the surface potential as 


$sheii('y,eo,a) 


Airaa i 

So 


E 

n^even 


2n -|- 1 
2 


X 


Pn 




In{eo)Pn{cOSv). 

(B14) 


Appendix C: Coulomb energy of a uniformly-charged disc 

In this appendix, we derive the exact expression for the 
Coulomb energy of a uniformly-charged circular disc of ra¬ 
dius a, total charge Q, and uniform charge density a = 
Q/{'Ka?'). The potential on the surface of the disc as a func¬ 
tion of p, the radial coordinate, has been derived in Ref. 13 
and is given by: 


V{p) = 4aaE , (Cl) 

where E{m) is the complete elliptic integral of the second 
kind: 



(C2) 


Note that 0 < P< a. This result can be used to obtain the 
electrostatic energy [/disc of this disc. Interestingly, while the 
potential on the disc surface is only available as an elliptic 
integral, we will soon see that the total electrostatic energy of 
the disc reduces to a simple form. 

Starting with the definition of the electrostatic energy, U = 
(1/2) J aVdA, we have 


Ddisc = / V{p)2-Kpdp. 

^ JO 


(C3) 


Substituting V{p) from Eq. (Cl) and using Eq. (C2), we ob¬ 
tain 




(/disc = dTTCT a / dp 

Jo 


d9\ll — ^ sin^O p. 


(C4) 


Carrying out the integral with respect to p first by employing 
the substitution t = \ — p^Av^djaJ, we obtain 


„3 

[/disc = dTTij"— / 

b Jo 


1 — COS^0 

sin^0 


dQ. 


(C5) 


The integral over 9 equates to 2 leading to the following ex¬ 
pression for the energy: 


The shell surface area element in the oblate spheroidal co¬ 
ordinates is given by 


TT StT 2 3 

(/disc — ^ ' 


(C 6 ) 


dA = sinu dudi/. 


To our best knowledge, this expression is reported here for the 
(B15) first time. 
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